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ABSTRACT 


Problems associated with subsynchronous motion at approximately 450 Hz for 
FPL turbopumpe Is the subject of this study. The basic model for the HPOTP 
is discussed including the structural dynamic model for the rotor and housing, 
component models for the liquid and gas seals, turbine-clearance excitation 
forces, and impeller-diffuser forces. Linear and nonlinear analyses are 
carried out to meet the following objectives: 

(b) Provide an explanation for the observed subsynchronous motion via 
nonlinear simulation results. 

(c) Evaluate proposed hardware changes in the HPOTP to remedy current 
difficulties . 

The analysis results support the following conclusions: 

(a) The current HPOTP is marginally stable and is subject to excessive 
synchronous bearing loads due to its proximity to the second 
critical speed. 

(b) Observed subsynchronous motion can be s'mulated by transient models 
which include bearing-clearance nonlinearities. Tl.^se nonlinearities 
generate forced subharmonic response at observed frequencies for the 
models considered. Bearing clearances and very light damping are 
sufficient to develop subsynchronous motion in an otherwise (zero - 
bearing-clearance) stable model. 

(c) Replacing the boost-impeller labyrinth seals with "damper" seals 
without stiffening the rotor reduces synchronous bearing loads, and 
Improves stability, but it does not necessarily eliminate subsynchronous 
motion. 

(d) By linear analysis, replacing the current main impeller with a shrouded- 
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inducer impeller eliminates both synchronous and subsynchronous 
problems associated with the second critical speed. 

(e) Replacing the current rotor with a stiffened rotor without 
changing seals eliminates subsynchronous motion but does not 
necessarily reduce synchronous bearing loads. 

(f) Replacing the current rotor with a stiffened rotor and replacing 
the current boost impeller labyrinth seals with damper seals 
sharply reduces synchronous bearing loads, Increases rotor stability, 
and eliminates subsynchronous motion. 
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I . INTRODUCTION 


A discussion of the rotordyr .unlc characteristics of the HPOTP is expedited 
by a review of figures 1, 2, and 3 which illustrate, respectively, the SSME 
powerhead and turbopump arrangement, the rotating assembly, and two of the 
zero-running-speed, coupled, rotor-housing modes fo.- the current (unstlf fened- 
rotor) turbopump. The two modes of figure 3 involve relatively small motion 
of the housing with large rotor motion corresponding to the first and second 
bending modes of the rotor. The 179 Hz mode is dominated by the overhung- 
turblne displacement, while the maximum displacement of the 530 Hz mode is 
at the main impeller, with relatively small motion at the turbines. The 
modes of figure 3 are calculated for the bearing stiffnesses = 500,000 
lbs/in. The calculated linear critical speeds for the turbcpump including 
gyroscopic coupling, seal forces, Impeller forces, etc. are at approximately 
12,500 cpm (208 Hz) and 32,500 cpm (540 Hz). 

Initial rotordynamlcs problems with the HPOTP involved subsynchronous 
motion at approximately 200 Hz associated with the first critical speed. This 
problem was eliminated by changing the turbine interstage seal from a stepped- 
labyrinth configuration to a constant-clearance, honeycomb-stator seal. The 
turbopump still displays subsynchronous motion in the 250 Hz frequency range 
at start up, and subsynchronous motion associated with the first critical 
speed continues to be of concern at higher power levels. 


The second problem associated with the HPOTP involves subsynchronous 
vibrations at approximately 400-420 Hz. This motion tended to initiate when 
the running speed reached the appropriate excitation frequency, and persisted 
as subsynchronous motion when the running speed was increased. A failure and 
explosion of the HPOTP occurred in a test run when the running speed first 


V 


“"SSw Q® 



Figure 1. SSME powerhead component arrangement and local 
rotor coordinate systems, 
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Figure 3. Undamped, zero-running-speed , nominal-rotor/housing 
modes '.ssociated with the first and second critical 
speed. 
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traversed the maxlmum-ampiltude frequency range and then was reduced to and 
remained at the now synchronous frequency, The failure was attributed to 
an inadequate bearing-carrier design which caused unequal load sharing of 
the bearings. Subsequently, redesigned bearing carriers eliminated the 
"400 Hz" subsynchronous motion. The general feeling w.s that the new 
bearing carrier designs caused n more equal load sharing and maintained 
the correct axial bearing preloads, thus realizing the "expected" radial 
stiffness of the bearings and elevating the rotor second critical speed. 

The bearing-carrier redesign was sufficient to yield a reliable pump 
for RPL operations; however, demonstration of FPL capability has been impeded 
by repeated occurrence of subsyr-^hronous motion in the 450-500 Hz frequency 
range. The incidence of this motion rarely occurs during initial operation 
at higher power levels, but tends to develop after an accumulation of time 
and develops more rapidly with higher synchronous vibration levels. However, 
once a unit begins to "whirl", it will repeatedly whirl, generally with 
progressively increasing severity. The frequency of the subsynchronous 
motion is sometimes at a fixed fraction of running speed ranging from 88 to 
95%. This "tracking" characteristic is in marked contrast to the earlier 
"400 Hz" phenomenon where the whirl-frequency of the subsynchronous motion 
remained constant, relatively independent of running speed, 

Post-run inspection of turbopumps has revealed that, although the balls 
in some bearings, primarily bearing 2, appear to be unblemished, their 
diameters have in fact, been reduced by as much as 0.17 nun in some cases. The 
cause of this dimensional attrition is unknown, but has been ascribed to either 
excessive radial loads or the absence of radial loads which leads to skidding. 

At various times, the following hypotheses have been put forward to 
explain the observed "450 Hz" subsynchronous motion: 


t 


i 


(a) Motion associated with the second-critlcal-speed modeshape of 
figure 2 is marginally stable, with the destabilizing forces 
provided by impeller-diffuser interaction forces [1], The 
principal deficiency of the linear-stability-analysis approach 
is the predicted whirl frequency of 530 Hz versus observed 
whirl frequencies around 450 Hz; hence, a linear model with 
nominal bearing coefficients predicts "supersynchronous" motion 
Instead of subsynchronous motion. The bearing stiffness coefficients 
of linear models must be reduced substantially below current 
estimates to yield subsynchronous motion at the correct whirl 
frequency. 

(b) The subsynchronous motion arises from bearing-stiffness asymmetry 
which rotates at the cage precessional speed and yields an 
additional forced excitation at approximately 90% of running speed. 
This is the only hypotheses which explains the "tracking" observed 
in the data and was put forward by Bernle Rowan at Rocketdyne 
before his untimely death. The problem with this explanation is 
that, according to Rocketdyne, backward-whirling motion would be 
excited, and the backwards modes are predicted to be heavily damped, 
The phenomenon of bearing-stiffness asymmetry exciting subsynchronous 
motion had been investigated by Yamamoto [2], who demonstrated 
experimentally that insertion of a single large ball into a ball 
bearing could yield stable, forced-subsynchronous motion, 

(c) The subsynchronous motion is a nonlinear phenomenon arising from 
motion in and out of the bearing dead-band* clearances. To the 

*Clearances provided between the outer case of the bearing and the housing to 
allow axial slipping. 
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author's knowledge, Steve Winder at MSFC first suggested this 
mechanism following the "400 Hz" HFOTP explosion. Simulations 
by the author [1] show that the bearing dead-bands can lower 
the peak vibration amplitudes into the operating range, as 
predicted by Yamamoto [3]. However, no demonstration has been made 
that the bearing clearances yield a lower whirl frequency than predicted 
by linear analysis. 

The present report aocuments the results of (a) linear analysis to characterize 
the nominal model and (b) nonlinear transient simulations which support the 
hypotheses that bearing dead-band clearances are principally responsible for 
the observed subsynchronous motion . No analyses were carried out to examine 

Rowan's hypothesis of rotating assymetry. 

While the causes of subsynchronous motion remain uncertain, the following 

har' • • modifications are being undertaken and considered: 

(a) The rotating assembly is being stiffened in the vicinity of bearings 
1 and 2. 

(b) The current labyrinth designs for the boost-impeller wear-ring 3eals 
are being replaced with damper seals which use a roughened stator and 
a smooth rotor. A swirl brake to reduce the inlet tangential velocity 
is also being implemented for the inlet seal. 

(c) Replacement of the current unshrouded inducers with shrouded inducers 
has been suggested, with the seals at the outer surface of the inducer 
developing relatively large stiffness and damping. 

The consequences of these changes on the rotordynamlc characteristics of the 
HPOTP are examined in this study . 

Continuing uncertainties exist as to the correctness of models for the 
Individual force elements of the rotordynamlc model. The author carried out 
a systematic examination of the literature concerning "best" component models 
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at the start of this study [4], and the preceding Phase B report [5] utilised 
the best available models and best available current estimates for parameters 


of these models. The current report provides a review of the changes and 


advances which have been made in component models of seal forces and again 



In summary, t h e results reported in this study were carried out to 
meet the following objectives : 

(a) develop the best-available model for the HPOTP in terms of forces 
developed by seals, impellers, turbines,' etc., 

(b) examine the linear dynamic characteristics of the HPOTP in terms 
of stability and synchronous-response predictions, 

(c) evaluate and develop explanations for the observed "450 Hz" sub- 
synchronous motion, and 

(d) evaluate the prospects of proposed hardware modifications for 
improving the dynamic characteristics of the HPOTP. 

The following chapter briefly describes the structure and components of the 
rotordynamlc model, with succeeding sections describing the results of 
linear and nonlinear analyses. 
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II. THE ROTORDYNAMICS MODEL 


Introduction 

This section briefly reviews the elements of the HPOTP rotordynamic models 
used in this study. Some of the actual "numbers" used to define the model 
are provided In Appendix A. 

Structural Dynamics Model 

The structural dynamics model of the rotor and the housing form the basic 
framework of the turbopump model, Both the housing and rotor structural 
dynamics models used here were developed by Rocketdyne. A general three- 
dimensional finite-element approach was used to model the IIPOTP housing 
without the rotor. A lumped parameter model was developed for the rotor 
using beam structural elements and lumped masses and inertias, The rotor- 
dynamics model uses modes from the housing model (without the rotor) and 
free-free rotor modes. The current (unstiffened) rotor model was developed 
by Bernie Rowan several years ago. The stiff ened-rotor model was provided 
in 1983 by Robert Beatty. 

Bearings 

The bearings are the structural elements which tie the rotors and housings 
together. As illustrated in figure 2, there are two sets of bearings. The 
net load from each set is transmitted through a plate-cylinder structure to 
the housing. The bearings in a bearing set are axially preloaded against 
each other, but are not designed to accept axial thrust loads from the turbo- 
pumps. A balance-piston arrangement at the discharge of the main impeller 
absorbs axial thrust, and radial clearances are provided at the bearing outer 
races to allow the bearing? to slip axially without developing excessive 
axial loads. 
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Experience and limited test data [6] for the SSME bearings indicate 
a nominal stiffness of approximately .5 X 10° lbs/in which is approximately 
one half of the values predicted by A. B. Jones-based analysis [7]. This 
nominal value was used throughout the current study. Prior analyses [5] 
examined the influence of feasible changes in the bearing stiffnesses and 
show the expected results, viz., the first and second critical speeds move 
up and down as the bearing stiffnesses are increased or decreased. 

The bearings are in series with the intermediate structure which 
transmits their load to the housing. The values used for those local 
stiffnesses are 

K , = K 0 = 2 X 10 6 lbs/in 

sl s2 a) 

K - ■ K , = 4 X 10 6 lbs/in 

S J S 4 

The radial clearances provided at the bearings to permit axial motion 
provide an essential "dead-band" nonlinearity. The clearance values used 
throughout most of this study are 

6 p = 6 = .5 X 1(T 3 in (2) 

where P and T denote pump and turbine bearings. 

Liquid Seals 

Liquid wear-ring seals are provided at the inlet and discharge of the boost 
impeller. The current inlet seal is a stepped labyrinth design with four 
cavities. The current discharge seal is a three segment, stepped seal. Each 
constant-radius seal segment has a series of circumferential grooves. 

Experience and limited test data [8] have shown that labyrinth or serrated 
seals of the type currently employed on the impeller inlet and discharge yield 
stiffness and damping coefficients which are substantially smaller than 
corresponding values for smooth constant clearance seals. Replacing the current 
grooved boost-impeller wear-ring seals with plain a.i. ,r seals has the potential 
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for a beneficial increase in stiffness and damping in the 11P0TP. In fact, 
as the following discussion explains, various additional possibilities exist 
for optimizing seals. 

The force-motion model for liquid seals has the form 
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In this model (X,Y) are the components of the relative motion between the rotor 
and housing. Black, et al. [9, 10] were responsible for most of the analytical 
developments related to the analysis of seals leading to the definition of 
stiffness, damping, and added-mass coefficients. His analysis demonstrates 
that the "cross-coupled" stiffness coefficient, k, arises sole’y due to fluid 
rotation within the seal. As a fluid element proceeds axially along an 
annular seal, shear forces at the rotor accelerate or decelerate the fluid 
tangentially until an asymptotic value is reached. For a seal with the same 
dlrectionally-homogerieous surface-roughness treatment on the rotor and the 
housing, the average asymptotic tangential velocity is Raj/2 where R is the 
seal radius and U) is the rotor running speed. 

The cross-coupled stiffness coefficient k acts in opposition to the 
direct damping coefficient G to destabilize rotors. Hence, steps which can 
be taken to reduce the net fluid rotation within a seal will improve rotor 
stability by reducing k. The modifications which have been undertaken for 
the boost-pump impeller seals propose to increase net damping by the following 
physical mechanisms: 

(a) The analysis of von Pragenau [11] has recently predicted that the 
asymptotic tangential velocity can be modified if a different 
surface roughness is used for the rotor and stator elements. Rough 
rotor/smooth stator and smooth mtor/rough stator combinations yield 
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higher and lower asymptotic values, respectively. Von Pragenau 
calls the rough-stator/smooth-rotor configuration a "damper seal" 
because of its enhanced stability characteristics. The rouglness 
has the additional potential benefit of reducing leakage. Damper 
seal configurations are presently under consideration for both the 
inlet and discharge seals of the boost impeller. Analysis by Childs 
and Kim [12] based on a finite-length solution [13] yield predictions 
which are consistent with von Pragenau. Rocketdyne develops the 
desired roughness on both boost-impeller seals by a knurled 
indentation pattern on the stator. 

(b) If the inlet tangential velocity can be reduced, the high axial 

velocities in a cryogenic seal are such that the fluid may proceed 
through a seal without substantially increasing its tangential 
velocity. An anti-vortex web has been introduced at the inlet to 
the boost-impeller inlet seal to reduce the inlet tangential 
velocity and yield a reduced k. This practice has been followed 
previously for the labyrinth seals of high pressure compressors [14] . 

Tests of an annular seals in which the stator has the knurled roughness 
pattern used by Rocketdyne and a smooth rotor have recently been carried out 
at TAMU. The results, which are documented in [15], show good agreement 
between theory [12] and experiment with respect to k and G, but a serious 
underprediction of K and M. For the present study, calculated values are 
used for k, C, m, and M, but calculated values for K are corrected upwards 
based on the results of [15]. For calculation, inlet tangential velocity for 
the discharge seal is assumed to be asymptotic, i.e., Rw/2. However, the 
swirl web at the inlet seal is assumed to be effective in reducing this value 
by a factor of two. Seal coefficients used in the rotordynamlc calculations 
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are provided in Appendix A. 

Gas Seals 

The HPOTP turbines are shrouded, and single-cavity tip seals are provided to 
reduce leakage between the turbine shroud and the stator. The interstage 
seal between the turbines uses a honeycomb stator element with a smooth rotor 
and inlet anti-vortex ribs to reduce the inlet tangential velocity component. 

A floating-ring shaft seal is provided at the turbine discharge to restrict 
leakage of the hot turbine gases towards the lox within the pump. The flow 
across this seal is choked. 

Rotordynamic coefficients for the tip seals use Rocketdyne calculations. 
Nelson's [16] method is used for the turbine interstage seal and the high- 
pressure turbine seal. His analyses apply to constant-clearance or convergent 
taper geometries, account for the development of tangential velocity within 
the seals, and different but direcfionally-homogeneous surface roughness on 
the rotor and stator. 

A constant-clearance configuration was originally used for the turbine 
interstage seal. More recently, Rocketdyne has replaced the constant- 
clearance configuration with a convergent taper. Computed rotordynamic 
coefficients are provided in Appendix A for both configurations, 

The "floating-ring" nature of the high-pressure turbine seal is 
such that its full reaction force can not be transmitted to the pump housing. 
Specifically, the seal is designed to slip before any appreciable force can be 
transmitted. To compensate for this situation, the calculated seal rotor- 
dynamic coefficients are arbitrarily reduced by a factor of 5. Calculated 
values are provided in Appendix A. 
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Turbine Clearance Excitation Forces 


Clearance-excitation forces are developed by turbines due to the dependency 
of local efficiency on local clearances. The destabilizing force is modeled 
by 


F x, 

1 ■ 

'0 k T " 


IX 
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(4) 


where T is the turbine torque, is the average pitch diameter of the turbine 
blades, H is the average height, and 3 defines the change in turbine efficiency 
due to uniform changes in clearance. Again, the components (X,Y) of Eq. (4) 
define the displacement of the turbine relative to the housing. Thomas initially 
identified this destabilizing phenomenon [17], while Alford [18] subsequently 
and independently developed the same model. Test results for shrouded turbines 
have yielded values for 8 on the order of 0.6 [19]. The value of 0.4 is used in 
this study. 


The dimensions of the JIPOTP turbines are 

D = 9.57 in, H = .496 in 
P 

The torque and clearance-excitation coefficient (8 = 1) are listed below 



co (cpm) 

T (in -lb) X 10 4 

1^ (lb /in) 

FPL 

30,960 

6.20 

13,060 

MPL 

19,841 

2.48 

5,230 

Note that the K T values 

are for 

both turbine stages with 

3 = i. 
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Impeller-Diffuser Forces 

A test program has been under way at the California Institute of Technology 
for some time to measure the static and dynamic forces experienced by a 
pump impeller in either a volute or a vaned diffuser, Chamieh et al. [10] 
have defined the following model on the basis of static measurements of an 
impeller within a vaned diffuser: 
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where is the Impeller radius, p is the fluid density, V2 = ^2 U is the impeller 
tip velocity, and A2=2iTR2b2 is the exit flow area. Note that the direct- 
stiffness coefficient in Eq. (5) is negative, i.e., the impeller-diffuser 


force causes a reduction in rr tor stiffness 
impeller-diffuser coefficients are defined by 

P A 2 V 2 

K = K* ~ 2 R - = K* (TT P b 2 R 2 2 )0J 2 

k = k*(TTpb 2 R2 2 )w 2 


From Eq. (5), the dimensional 


( 6 ) 


More recently. Cal. Tech, researchers [21] have reported preliminary results 
for an extension of the model provided by Eq. (5) to include direct and cross- 
coupled damping terms. However, rotordynamic analysis of the HPOTP suggests 
that the proposed values are unrealistically high. In particular, the 
proposed direct damping value is large enough to eliminate any possible 
rotordynamic problems. The cross-coupled damping coefficient acts like a 
"gyroscopic-stiffening" element an’, also seems to be too large, since it 
would elevate the predicted whirl frequency to 638 Hz. In view of these 
uncertainties, the original simple model of Eq. (6) was used for rotor- 
dynamic calculations with the damping model of [21] discarded. 
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The dimensions of the two impellers are provided below: 


Main Impeller: 1*2 = 

3.35 in, b 2 = 

1 in 


Boost Impeller: « 

2.60 in. b 2 a 

.270 in 


The following coefficients result 




Main Imneller 

Boost Impeller 

FPL 

MPL 

FPL 

MPL 

w (cpm) 30,960 

19,841 

30,960 

19,841 

P (Kg/m 3 ) 1,137 

1,137 

1,114 

1,114 

K (lb/in) -7.88x10" 

-3.24X10'* 

-1.28X10" 

“5 , 2‘ 

k (lb/in) 2.76X10'* 

1.13X10" 

4,480 

1,840 

Fixed-Direction Side Loads 




The hydrodynamic side loads are assumed to be proportional to speed 

The proportionality constants are 

listed below 




K (X - Z) 

K (Y 

- z) 


lbs • sec 2 

lbs • 

sec 2 

Boost Impeller 

-2.806 X 10“ 

5 4,950 

X 10" G 

Main Impeller 

7.281 X 10“ 

5 6.117 

X 10“ 5 


Turbines 


2.379 X 10 


“5 


These coefficients yield a total main-impeller side load of 960 lbs at FPL. 
The side load is used in the nonlinear analysis and, by its interaction 
with the dead-band clearance and Imbalance magnitude, has a significant 
influence on rotor response. 
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Balance-Piston Stiffness and Damping Coefficients 


The balance-piston stiffness is modeled by the quadratic relationship 
K„ ° -431.97 u + .3542 w 2 (7) 

where w is the running speed in rad/aec. However, K z is never allowed 
to fall below 200,000 lbs/ln. Eq. (7) fits Winder's graphical data [22], 
Balance-piston damping is held at 15% of critical for all speeds. Axial 
motion of the rotor is only coupled to the housing in the nonlinear model. 
The linear model does not include this feature. 

Imbalance Distribution 

The imbalance distribution used in all cases consisted of the following 


aligned imbalances: 

Locaclon 

Magnitude 

(a) 

Boost Impeller 

.1273 gm cm 

(b) 

Main Impeller 

10.18 gm cm 

(c) 

Mid- turbine 

12.73 gm cm 

Table 

1. HPOTP Imbalance 

distribution. 


While considerable uncertainty exists concerning the particular Imbalance 
distribution in a given turbopump, the distribution of table 1 provides 
adequate excitation for the modes of interest. 
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LINEAR ANALYSIS 


Introduction 

As noted above, bearing "oead-band" clearances are the significant known 
nonlinearity which Influences roto.dynamlcs . If the dead-band clearances 
are neglected, a linear model results. The present chapter concerns 
results of '’•'.near analysis for various models of HPOTP turbopump configurations. 
The following chapter concerns nonlinear-analysis results which differ in 
many respects from linear predictions. The linear analysis results of the 
present section provide an efficient and helpful characterization of the 
turbopump’s rotordynamic characteristics. 

The analysis procedure used here is bosically the same as that outlined 
in reference [23]. Modal coordinates based on the zero-running-speed coupled 
rotor-housing modes are used. Gyroscopic coupling and forces due to seals, 
turbine clearance-excitation, the interaction of impellers and diffusers, 
damping, etc. couple the modal coordinates via modal stiffness, damping, and 
inertiu matrices. The onset speed of instability for a turbopump configuration 
is defined by calculating the complex eigenvalues of the system dynamic matrix 
at various speeds. Synchronous-response amplitudes of bearing reactions and 
acceleration levels of accelerometers mcuuted on the turbopump housing due to 
imbalance are calculated. 

Table 2 provides the distinguishing characteristics of the six turbopump 
configurations which are analyzes in this section. Except as otherwise noted, 
the six turbopump models are identical and use the nominal parameters of the 
preceding section and Appendix A. 

Parameters of the "nominal" model were selected to represent the current 
flight configurations, 'without additional damping, the nominal model predicts 



Rotor 

Boost-Impeller 
Wear-Ring Seals 

Main 

Impeller 

a. Nominal Model 

current 

Labyrinth 

Unshrouded 

b. Nominal Model 

with Damper Seals 

current 

Rocketdyne 

Damper 

Unshrouded 

c. Nominal Model 

with Shrouded Inducer 

current 

Laybrinth 

Shrouded 

d. Stiffened Rotor 

stiffened 

Labyrinth 

Unshrouded 

e. Stiffened Rotor with 
Damper Seals 

stiffened 

Rocketdyne 

Damper 

Unshrouded 

f. Stiffened Rotor with 
Otimized Damper Seals 

stiffened 

"Optimized" 

Damper 

Unshrouded 


Table 2. Turbopump configurations for linear analysis. 
























an OSI sunset Speed of Instability) at 17,200 cpra. Hence, the model was 
"tuned" by adding concentrated linear damping at the main Impeller. This 
damping develops forces which are proportional to the relative velocity 
between the rotor and housing at this location. The linear damping 
coefficient, C a 5 lb sec/in, was sufficient to elevate the OSI to 30,480 cpm 
and is used for all models. No damping is provided at the bearings. The 
configurations of table 2 were selected to characterize the current flight 
configurations and examine the consequences of the following changes to this 
configuration: 

(a) Replace the labyrinth, boost-impeller, wear-ring seals with damper 
seals . 

(b) Replace the unshrouded Inducer with a shrouded Inducer. 

(c) Stiffen the current rotor. 

(d) Stiffen the current rotor and use damper seals for the boost-impeller 
wear-ring seals. 

(e) Stiffen the current rotor and use "optimized" damper seals for the 
boost impeller. 

Results for these configurations are provided below. 

Current Flight Configuration 

Figure 4 Illustrates the bearing reactions for running speeds from 5,000 to 
40,000 cpm. The first and second critical speeds are clearly evident, with 
the second critical speed having very little damping. Figure 5 (a) illustrates 
the same results for speeds from 5,000 cpm to FPL. Figures 5 (b) and (c) 
illustrate predicted acceleration levels for accelerometers mounted on the 
housing at the pump and turbine ends of the turbopump. Clearly, FPL is quite 
near the predicted second critical speed. Also, while the bearing reactions 
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SEARING I 
SEARING 2 
BEARING 3 
SEARING 4 


Figure 4, Bearing reactions for the nominal (current) 
linear model with speeds out to 40,000 cpm. 
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BEARING I 
BEARING 2 
BERRING 3 
BEARING 4 


Figure 5(a). Bearing reactions for the nominal linear model 
with speeds out to FPL, 
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ORIGINAL RAGS fg 
OF POOR QUALITY 


FREQUENCY RESPONSE ANALYSIS SSME HPOTP 
NOMINAL BASELINE MOOEL 

W 1 0 • T 



RUNNING SPEED (RPMI 


X*OI RECTI ON 
Y'DIRECTION 


Figure 5(b). Pump accelerometer levels in the X-Z and Y-Z 
planes . 


FREQUENCY RESPONSE ANALYSIS SSME HPOTP 


NOMINAL BASELINE MODEL 



RUNNING SPEED (RPMI 


X-OIRECTION 


Y-DIRECTION 


Figure 5(c). Turbine accelerometer levels in the X-Z 
and Y-Z planes. 
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are responsive to the "rotor" critical speeds, the housing accelerometer levels 
may respond more sharply to rotor-housing combined modes. This is particularly 
true for the first critical speed at approximately 12,500 cpm. However, the 
predicted accelerometer levels are quite responsive to motion associated with 
the "rotor" second critical speed at approximately 32,500 cpm. 

The above results support the following conclusions with respect to the 
dynamic characteristics of the HPOTP: 

(a) Response . The proximity of the second critical speed to FPL is a 
continuing cause for concern with respect to excessive bearing 
loads. Even modest losses of bearing stlf fn^canq are sufficient to 
drop the second critical speed into the operating range. 

(b) Stability . Based on the prior "400 Hz" experience with this turbopump, 
the second critical speed appears to be lightly damped and subject 

to instability. 

Changes in the turbine-interstage seal coefficients for the current 
(convergent-tapered) and prior (constant-clearance) configurations made no 
appreciable, difference in the rotordynamic characteristics of the nominal 
model with respect to second-critical speed response. 

Seal Modifications with the Current Rotor 
BooAt ImpeZtoA WnoA-Zing Szati, 

Figure 6 illustrates the predicted bearing reactions which results with damper 
seal configurations for the boost impeller, and a comparison of this result 
with the nominal-rotor-model results of figure 5 (a) indicates a sharp reduction 
in peak bearing-reaction magnitudes. Additionally, the OSI is increased from 
30,480 to 36,350 cpm. 

Sh/Lc uded I nduceJL Seciti 

Figure 7 illustrates the results for the shrouded inducer configuration. Stated 


BERRING RFRCTICNS 6ERRING RERCTIQNS 


ORIQWtt Tt- v : 

OF ROOK ; . ’■/ 


FREQUENCY RESPONSE 6SME HPOTP 
FWR & RWR DRMPER SERLS 



RUNNING SPEED "RPMI 


BEARING I 
BERRING 2 
BERRING 5 
BERRING 4 


igure 6, 


Bearing reactions for the nominal linear 
model with "damper" seal coefficients for 
boost-impeller seals. 


FREQUENCY RESPONSE ANALYSIS SSME HPQTP 
SHROUDED INDUCER SERLS 



RUNNING SPEED 1RPM) 


BERRING I 
•BERRING 2 
BERRING 3 
BERRING 4 


Figure 7. Bearing reactions for the nominal rotor 

model with "dan, per" seal coefficients for 
a shrouded inducer. 


- 22 - 



briefly, the second critical speed Is simply eliminated i'or this configuration, 
and the OSI is Increased to 59,950 cpm. 

Stiffened Rotor Results 

Figure 8 illustrates rotor-housing mode shapes which are dominated by the rotor 
motion and correspond to the first and second rotor critical speeds. Observe 
that the modeshapes are similar to those of figure 3 but that the natural 
frequency associated with the second critical speed has been increased from 
530 Hz to 558 Hz, This increase is the expected result associated with an 
increase in rotor stiffness. However, the corresponding result for the over- 
hung turbine mode is a decrease from 179 Hz to 168 Hz, which is certainly not 
expected. Bob Beatty at Rocketdyne feels that the results for the original 
rotor model are in error because of a low mass estimate in the turbine disks. 
Specifically, the mass due to goldplating was not accounted for in the original 
analysis. 

StL^&ne.d-Rotoa/Baoit-lmpaLteA.-Labtj^uvth-Se.ccl Model 

Figure 9 illustrates the bearing-reaction-magnitudes versus speed results when 
the nominal rotor structural-dynamic-model data is replaced with the stiff-rotor- 
model data. Observe that the second critical speed is elevated from approximately 
32,500 cpm to 35,000 cpm, and that the FPL bearing loads are reduced. The OSI 
is increased to approximately 40,000 cpm. 

Sti^e.n&d-Rotoa/Boa6t-lmpQ,tieA. Vamp aji Seal 

Figure 10 illustrates the bearing-reaction-magnitude versus speed results for 
a stif fened-rotor/damper-seal configuration. By comparison to figure 9, the 
second critical speed is elevated to 42,000 cpm, and the peak bearing-reaction 
magnitudes are decreased. For this configuration, the OSI is increased to 
approximately 45,000 rpm. 

Sti^m&d RotoA/sie.duc.e.d-StL^n&AA Boo4t~7mpe££eA VampeA Se.at6 

The predicted direct stiffness of a damper-seal configuration for the boost- 
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Figure 9. 


Bearing reaction magnitudes versus running speed for the 
stiff ened-rotor model with labyrinth boost-impeller seals 
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Figure 10. Bearing-reaction magnitudes versus running speed for the 
stif fenod-rotor model with damper seals for the boost 
impeller. 

HPOTP STIFF ROTOR FREQUENCY RESPONSE 
say. KXX FWR » RWR damper SEALS 



5000. 15000. 25000. 35000. 

10000. 20000. 30000. 40000. 
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BEARING 1 
BEARING 2 
BEARING 3 
BEARING 4 


Figure 11. Bearing-reaction magnitudes versus running speed for the 
stif f ened-rotor model with reduced-stiffness damper seals- 
for the boost impeller. 
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impeller discharge seal is comparable to that of the bearings and acts to 
markedly reduce the amplitude of relative rotor motion at the sea] location. 
Beyond a certain point, increasing the seal's direct stiffness becomes 
counterproductive in that the effectiveness of the seal's damping is limited 
in reducing synchronous-response amplitudes. Tests of seal stator configura- 
tions [15] suggest that seals can be developed with increased net damping and 
substantially decreased direct stiffness coefficients. The results of figure 
11 illustrate the consequences of a 50% reduction in the predicted direct 
stiffness coefficient of the boost-impeller damper seals. While the bearing 
reaction magnitudes are relatively unchanged. The OSI for the arrangement is 
sharply elevated to 62,000 cpm. 


SUMMARY AND CONCLUSIONS 

Table 3 contains a summary of the pertinent results for linear analysis. 
An assessment of the results presented in figures 6 through 11 and table 3 
supports the following conclusions: 

CuJUiant R o toh. 

(a) Modification of the boost impeller seals has the potential for 
significantly reducing bearing reactions at FPL; however, the 
second critical speed remains only slightly above FPL and any 
loss of bearing stiffness will drop it into the operating range. 

(b) The damper-seal modification yields a predicted increase of the 
OSI by 19%. 

(c) The shrouded-inducer design eliminates both the second critical 
speed and the stability problem. 
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Conf iguration 

2nd Critical 
speed (cpm) 

OSI 

{’»: pm) 

Whirl 
Freq. (Hz) 

FPL Bearing Reaction (lbs) 

1 

2 

3 

4 

Nominal 

32,500 

30,480 

536 

553 

635 

667 

548 

Nominal with Damper 
Boost-Impeller Seals 

34,000 

36,350 

576 

251 

349 

513 

411 

Nominal with Shroud- 
ed Inducer Seals 

critically 

damped 

60,000 

80 

80 

85 

56 

58 

Stiffened Rotor 

35,000 

40,154 

583 

305 

314 

229 

173 

Stiffened Rotor 
with Damper Seals 

42,000 

45,000 

701 

64 

82 

90 

54 

Stiffened Rotor 
with 50% stiff- 
ness Damper Seals 

62,000 

44,000 

700 

82 

99 

102 

66 


Table 3. Linear OSI, whirl frequency, 2nd critical-speed location, and FPL bearing 
reactions. 
















































Sti^Q.ne.d Roto a 




(a) Stiffening Che rotor, without changing the seals, yields a 
substantial elevation of both the OSI (32%) and the second 
critical speed (29%) , and a sharp reduction in FPL bearing 
loads. Assuming that the structural dynamic models accurately 
reflect the relative stiffness of the two rotors, stiffening 
the rotor provides a marked improvement in the rotordynamic 
characteristics of the turbopump. 

(b) Boost-impeller damper seals provide a significant additional 
improvement in the dynamic characteristics of the I1P0TP, primarily 
in reduction of FPL bearing loads. 

(c) Reducing the direct stiffness of the boost-impeller damper seals 
mainly acts to elevate the OSI with minimal increases in FPL bearing 
loads . 
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NONLINEAR ANALYSIS PROCEDURES AND RESULTS 


Introduc t ion 

As noted previously, the bearing "dead-band" clearances provide the essential 
nonlinearity in the HPOYP model. The bearing clearances interact with the 
effects of side loads and rotor imbalance to yield significantly different 
results for a nonlinear model than those predicted by linear models. 

Figure 12 illustrates the analytical model used by Yamamoto [3] tn 
investigate the influence o F bearing clearances on rotordynamlc response in 
the absence of side loads. For zero bearing clearances, the model of figure 12 

reduces to a simple Jeffcott model with viscous external damping. Figure 13 

illustrates the response characteristics for a progressive increase in the 
ratio of bearing clearance to imbalance eccentricity £ = e/a. The results are 
for a damping ratio of 2.5%, and Indicate that the speed location jt maximum 

bearing reactions is reduced by increasing Moreover, the drop in amplitude 

for speeds above the maximum bearing-load can be precipltious . The response 
characteristics of figure 13 can give rise to "jump" phenomena with the 
synchronous vibration level jumping either up or down for very small changes 
in running speed. Flight data for the HPOTP have demonstrated sudden step 
increases in accelerometer levels. 

From these results, one would anticipate that bearing clearances could 
easily drop the peak-bearing-load running speed .location associated with the 
second critical speed into the operating range. In fact, parametric studies 
of the bearing clearances clearly confirm this result [1]. They also confirm 
that peak-bearing-load speeds can be reduced below FPL for sufficiently large 
clearances. Hence for a given running speed, a nonlinear model with bearing 
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Figure 12. Yamamoto's model for synchronous response with 
bearing clearances. 



Figure 13. Synchronous-response characteristics for 
Yamamoto's model with increasing bearing 
clearances . 
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clearances can yield either oubstantially larger or smaller bearing reactions 
than the linear model. 

The nonlinear simulation results of! this section were developed for the 
following purposes: 

(a) develop an explanation for the observed 450 Hz whirl frequency, 
and 

(b) evaluate the effectiveness of proposed hardware remedies for 
improving the rotordynamlc characteristics of the HPOTP. 

Current Rotor Configuration 

The nonlinear model was verified by comparison to results from the linear model, 
e.g., the second-critical-speed location at 32,500 cpm was obtained for zero 
bearing clearances. Figure 14 illustrates the bearing 2 reaction magnitude 
and turbine acceleration signals at FPL for the nominal model with the bearing 
clearances of Eq. (2). Observe the "clipping" of the bearing reaction magnitudes 
in figure 14 (a), which results from motion through the bearing dead-band. 

Motion in and out of the dead-band results as a combination of static displacement 
due to side loads and elliptical orbits due to Imbalance. Bearing clipping 
generates a nearly periodic impulsive loading at the bearings that excites the 
beating motion which is evident in the accelerometer signal prediction of figure 
14 (b) . Spectrum analysis results of the accelerometer signals at FPL and two 
higher speeds are provider in figure 15 and reveal a subsynchronous signal at 
26,500 cpm (443 Hz). Both the synchronous and subsynchronous signals are 
observed to decrease as the running speed increases. Observe that the frequency 
of the subsynchronous motion is consistent with frequencies observed in test 
data as apposed to the linear predictions of 530 Hz. The linear model corres- 
ponding to the result of figure 14 is lightly damped but stable. The nonlinear 
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Figure 14 (a) . Bearing 2 reaction magnitude versus time at FPL 
for the nominal nonlinear model. 


TRANSIENT ANALYSIS SSME HPOTP 
STEADY STATE AT 30,381 RPM 



Figure 14(b). Turbine X-Z plane acceleration magnitude versus 
time for the nonminal nonlinear model. 
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model remains In a limit cycle motion with subsynchronous components until 
the clearance-excltation-force coefficients at the main Impeller are Increased 
to 250% of their nominal values, but then diverge exponentially. Based on 
these observations, the subsynchronous component of figure 14 is deemed to 
be the result of bearing nonllnearltles and not the result of an Instability. 
However, the nonlinear sabharmonlc motion can only be developed for light 
rotor damping, and the rotor is certainly unstable for zero bearing stiffnesses 
which is the case when motion is within the dead-bands . 

The subsynchronous motion of figures 14 and 15 arises at speeds above the 
maximum synchronous response speed predicted by Yamamoto's model in figure 13. 

This statement is supported by the results in figure 16 which shows the bearing 
2 reaction magnitude predictions at 29,450 cpm. Observe that the synchronous 
bearing-reaction magnitudes are significantly higher than the subsynchronous 
amplitudes of figure 14 (a). 

A second nonlinear model configuration which yielded a substantial 
subsynchronous vibration component was obtained by doubling both the Imbalance 
magnitudes and the destabilizing force coefficients at the main impeller, 
adding damping at the bearing (2 lb s/in) , and increasing the damping at the 
main impeller to 15 lb s/in. This configuration is predicted to be linearly 
stable at FPL. Spectral analysis results for this configuration at FPL and 
two higher speeds are presented in figure 17. In this case, the subsynchronous 
component increases with running speed, while the synchronous component decreases. 

Figures 18 through 20 illustrate results for the nominal nonlinear rotor mode 
with damper seals at the boost impeller. By comparison to figure 16, observe 
that the damper seals significantly reduces the RPL synchronous loads and the 
constant side loads on bearing 2. By comparison to figure 14 (b) , figure 19 
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Figure 15, Turbine X-Z plane acceleration spectr at FPL 
and higher speeds for the nominal nonlinear 
model , 
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Figure 16. Bearing 2 reaction magnitude at 29,500 cpm 
for die nominal nonlinear model. 
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Figure 17. Turbine X-Z plane acceleration spectra at FPL 
and higher speeds for nonlinear model number 2. 
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Figure 18. Bearing 2 reaction magnitude at 29,500 cpm for 
the nominal nonlinear model with damper seals 
at the boost impeller. 
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Figure 19. Bearing 2 reaction magnitude at speeds above FPL for 
the nominal nonlinear model with damper seals at the 
boost impeller. 
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Figure 20. Preburner X-Z accelerometer spectrum at FPL for 
the nominal nonlinear model with damper seals at 
the boost impeller. 
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demonstrates that the damper seals significantly alter the nature of the 
subsynchronous motion, but tend to increase the "clipping" of the bearing 
reactions. Presumably, the clipping motion is increased because of the 
direct stiffness of the damper seals which tends to counter the side loads 
and center the bearings within their clearances. The spectra of figure 20 
demonstrate that the subsynchronous motion remains significant for the 
nominal rotor model with boost-impeller damper seals. 

Stif fened-Rotor Results 

The results of the preceding section examine the consequences of a change 
in the current boost-impeller wear-ring seals from a labyrinth to a damper-seal 
design with no stiffening of the rotor. Conversely, the frames of figure 21 
illustrate the consequences of introducing the stiffened rotor without changing 
the seals. Observe that the FPL bearing reactions and acceleration levels have 
actually Increased as compared to the nominal model of figure 14 which contradicts 
the linear predictions. Introducing the stiffened rotor elevates the second 
critical speed and eliminates motion in and out of the bearing dead-bands, Pre- 
sumably, at some higher running speed, the dead-band motion and associated sub- 
synchronous motion would return. Fourier analysis of the acceleration signal of 
figure 21 (b) shows synchronous and twice synchronous components but no sub- 
synchronous . 

At present, the stif f ened-rotor and boost-impeller damper seals are 
a "package-deal" in that all modified units will incorporate both changes. 

Figure 2 illustrates predicted FPL results for this configuration. Observe 
by comparison to figures 14 and 21 that the bearing -2 reaction magnitude and 
the turbine -X accelerations are both reduced substantially. Fourier spectrum 
analysis of the acccl. signal of figure 22 (b) reveals synchronous and twice- 
synchronous components, but no subsynchronous. 
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Concerning the vibration characteristics at RPL, the results of figure 23 
illustrate a sharp reduction in the bearing -2 reaction magnitude as compared 
to figure 16. Fourier analysis of the turbine acceleration signal of figure 
23 (b) reveals no subsynchronous components. 

One precautionary note is emphasized with respect to the encouraging 
bearing-reaction predictions of figures 22 (a) and 23 (a) . The model is set 
up so that the boost-impeller seals and bearings are exactly centered. Hence, 
the seals act to unload the bearing alternating and constant loads. However, 
if the centers are radially displaced the seal stiffness would act to increase 
the steady-state bearing reaction magnitude. 
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Figure 21 (a) . 


Bearing 2 reaction magnitude versus time at FPL for the 
stiffened rotor with labyrinth, boost-impeller, wear-ring 
seals . 
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Figure 21 (b) . 


Turbine acceleration level versus time at FPL for the 
stif fened-rotor with labyrinth, boost-impeller wear-ring 

seals . 
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Figure 22 (a) . Bearing 2 reaction magnitude versus time at FPL for the 
stiffened rotor with damper-seal designs for the boost- 
impeller wear-ring seals. 
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Figure 22 
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for the boost 



impeller , 

wear-ring seals. 
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Figure 23 (a). Bearing 2 reaction magnitude versus time at FPL for the 

s t if fened-ro tor /boost-impeller damper seal configuration. 
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Figure 22 (b) . Turbine X acceleration level versus time at FPL for the 

stiffened rotor/boost-impeller damper seal confuguration . 

- 41 - 


CONCLUSIONS 


Linear and nonlinear analyses of the HPOTP support the following 
conclusions : 

(a) By linear predictions, the second critical speed of the current 
(unstlf fened-rotor) HPOTP is only slightly above the FPL running 
speed and is lightly damped. Even a modest loss of bearing 
stiffness is sufficient to drop the second critical speed into 
the operating range and yield a prediction of excessive bearing 
loads . 

(b) Based on prior test experience, the HPOTP appears to be lightly 
damped . 

(c) The subsynchronous motion which is evident In test results can 

be obtained in a nonlinear model due to bearing clearances. Motion 
in and out of the dead-band clearances excites subharmonic motion 
for a stable but lightly damped model. 

(d) A whirl frequency of subsynchronous motion at 440 Hz is predicted 
by the nonlinear model for the current rotor which is consistent 
with test results but at odds with a linear prediction of 530 Hz. 
This whirl frequency reduction results solely from the bearing 
clearances, and does not require any stiffness reduction in the 
model at the bearings or elsewhere. 

(e) Bearing clearances can drop the peak-vibration running speed ?nto 
■ the operating range yielding a prediction of excessive bearing 

loads at RPL. 

(f) Incorporation of damper seals into the boost impeller seals with 
the current rotor markedly reduces predicted synchronous bearing 
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loads, but Increases bearing-reaction clipping and subsynchronous 
motion at FPL. 

(g) Introducing a stiffened rotor without changing the boost-impeller 
seals elevates the rotor critical pneed and eliminates subsynchronous 
motion at FPL; however, predicted bearing-reaction magnitudes at 

FPL remain high. 

(h) A stiff ened-rotor/boost-impeller wear-ring damper seal configuration 
eliminates subsynchronous motion in the operating range and markedly 
reduces synchronous bearing reactions. 
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APPENDIX A: INPUT DATA 


Cun ant Rotor: Eigenvalues and Damping Factors 

The rotor eigenvalues anu eigenvectors used here are baaed on a model 
by B. Rowan. The free-free eigenvalues used are listed below. 


h 

C3 

0 

X 7 

S3 

3724 Hz 

X 2 

S3 

0 

X 8 

C3 

4389 Hz 

X 3 

S3 

426 Hz 

X 9 

« 

6600 Hz 

X 4 

C3 

970 Hz 

X 10 

S3 

7397 Hz 

X 5 

S3 

1561 Hz 

X 11 

S3 

10396 Hz 

X 6 

S3 

2698 Hz 

X 12 

E2 

11916 Hz 


One-half percent of critical damping was used for modes three through twelve. 
Zero damping was used for modes 1 and 2. 

Stiffened Rotor: Eigenvalues and Damping Factors 

The rotor eigenvalues and eigenvectors used here are based on a 1983 
model provided by Robert Beatty of Rocketdyne. 


X 1 

= 

0 

X 7 

= 

3211 

Hz 

X 2 

= 

0 

X 8 

= 

3787 

Hz 

X 3 

SJ 

451 Hz 

X 9 

S3 

4722 

Hz 

X 4 

= 

1053 Hz 

X 10 

= 

5571 

Hz 

X 5 

= 

1779 Hz 

X 11 

= 

6879 

Hz 

X 6 

= 

2998 Hz 

X 12 

= 

7604 

Hz 


One-half percent of critical damping was used for modes three through twelve. 
Zero damping was used for modes 1 and 2. 

Housing Eigenvalues and Damping Factors 

The case eigenvalues and eigenvectors are based on a 1982 Rocketdyne 
structural-dynamic model. The eigenvalues used in the study are 
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cl 

G3 

45 Hz 

X r 
c6 

CJ 

351 

Hz 

c2 

□ 

86 Hz 

X c7 

a 

432 

11 7 

c3 

a 

111 Hz 

X c8 

a 

468 

Hz 

c4 

u 

300 Hz 

X 0 
c9 

o 

CO 

CO 

<r 

Hz 

c5 

C3 

310 Hz 

X cl0 

CJ 

542 

Hz 


One-half percent of critical damping was used for all housing modes. 

Seal Rotordynamic Coefficients 

Experience has shown that turbopump seals have a very significant 
Influence on rotordynamics . Because of the HPOTP modeshapes, the turbine 
seal” have a predominant Influence on the first-critical response and 
minimal influence on the second-critical response. Conversely, seals at the 
boost or main impeller have a significant Influence on motion associated with 
the second critical speed, but minimal Influence on motion associated with 
the first critical speed, To improve rotordynamic response, the following 
seal changes have been Implemented or considered: 

(a) The turbine Interstage seal was first changed from a stepped laby- 
rinth to a constant-clearance/honcycomb-stator configuration, and 
then to a convergent-taper /honeycomb-stator configuration. 

(b) The boost-impeller wear-ring seals have been changed (in test pumps) 
from labyrinth to constant-clearance seals with roughened stators. 

(c) Sealing surfaces are to be created at the outer surface of shrouded 
Inducers for the main impeller. 

Rocketdyne-calculated rotordynamic coefficient values (1977) are used for 
labyrinth seals, viz., the turbine tip seals and boost-impeller labyrinth 
seals. Rocketdyne-calculated values are also used for the shrouded- inducer 
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/■ 


seals. However, TAMU calculated values are used for the boost-impeller damper 
seals, the high-pressure-turbine, and the turbine interstage seals. The analysis 
of Childs and Kim [12] Is used for Incompressible-! luld seals and the analysis 
of Nelson [16] 13 used for gas seals. For completeness, Input data for these 
calculations are also provided. 

Tables A,1 and A. 2 contain the appropriate data for the boost-impeller sea_s, 
while tables A. 3 and A. 4 contain data for the high-pressure turbine seal and 
the turbine interstage seal. Data for the remaining seals are provided in 
table A. 5 and A. 6. 
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Soul 

Power 

D 

L 

C r 

u 3o 

Y 

M X iW 

AP 

w 

hovel 

(in) 

(in) 

(in) 


lb/ft 1 

lb nne/fc a 

pni 

cpm 

Inlot 

PPI. 

2.89 

0.5 

.006 

.25 

70.2 

.373 

2006 

30.367 

MPI. 

2.89 

0.5 

.007 

.25 

70.0 

.385 

970 

20,298 

Discharge 

FPL 

3.345 

0.87 

.004 

1.0 

70.2 

.373 

6505 

30,367 

MPL 

3.345 

0.87 

.0048 

1.0 

70.6 

.373 

3246 

20,298 


*Uq o a U 0o /(Rw/2) = Normalized inlet tangential velocity 
Turbulence Empirical Coefficients 


Stator: 

ms “ -0.136, 

ns - -.1357 

Rotor: 

mr a -0.25, 

nr ° ,079 


Table A.l Input data for boost-impeller damper seal calculations. 



Power 

K 

k 

C 

c 

M X 10 1 * 

m X 10 5 

w 

Seal 

Level 

lb /in 

lb/in 

lb soc/in 

lb aec/in 

lb eec a /in 

lb uec a /in 

lb/aec 


FPL 

1 

21,280 

30. 

1.28 

8.2 

- 1.44 

8.37 


KPL 

46,960 

(93,920) 

8,800 

18. 

0.72 

6.86 

- 2.16 

7.07 

Discharge 

FPL 

732,220 
(1.464 X 10 6 ) 

369,800 

280. 

14.8 

80.6 

20.00 

8.28 

HPL 



163. 

8.42 

67.3 

17.0 

7.57 


Table A. 2 Calculated rotordynamic coefficients and leakage for boost- 
impeller seals. 

Direct stiffness values in parenthesis were used in rotordynamic 
calculations and are twice calculated values. 



























ORIGINAL ;4 

OF POOR QOAUi * 



Turbine I. S. 

II. P. Turbino 

IfPL 

FPL 

MPL 

FPL 

w (cpra) 

20,283 

30,376 

20,283 

30,376 

P ln (pain) 

2,857 

5,540 

1,983 

3,499 

T i„ <•“> 

1,024 

1,498 

897 

111.9 

p ex (P nlfl) 

2,354 

4,261 

45.6 

63,4 

R Ideal gas 
8 constant 

532 

445 

593 

530 

(lbm/ft-aec) 

1.02 X 10" 5 

1.44 X 10" 5 

.92 X 10' 5 

1.11 X 10 -E 

apec. heot 
ratio 

1.385 

1.371 

1.4 

MMI 

R (in) 

2.B67 

2.867 

1.316 

■rasa 

L (in) 

0.995 

0.995 

0.318 

0,318 

U 6o (o)/Rw 

0.25 

0,25 

0.5 

mm 


Turbine I. S. Seal 

Constant (radial) clearance: C in ° C ox = .013 in 

Convergent taper: C^ n = 0,15 in, C ex = ,010 in 

Hirs* coefficient rotor: m “ -0,1634, n = .03757 

r r 

Hirs' coefficient atator: ra » -.002512, n - ,01534 

fl B 

H, P, Turbine Sea l 

Radial Clearance MPL: C in » .005, C qx » ,0016 

Radial Clearance FPL: C in *> .0048, C ex «■ .0013 

Hira 1 coefficient rotor: ro “ -.170, n a 0,040 

r r 

Hirs' coefficient atator: m = -0.143, n„ ■ 0.030 

a a 


Table A. 3 Input data for turbine gas seals. 



Power 

K X 10 -5 

k 

C 

c 

• 

w 

Seal 

bevel 

lb/ in 

lb /in 

lb /in 

lb / in 

i b /r p r 

Turbine H.P. 

FPL 

1.965 

4.140 

2.51 

.015 

.0335 

MPL 

1.016 

1.560 

1.39 

.005 


Turbine I.S. 

FPL 

2.02 

9.440 

11.3 

.067 

. i i ( * j n , 

.7427 

Const. Clearance 

MPL 

.918 

3.050 

6.56 

. 041 

.3000 

Turbine I.S. 

FPL 

1.05 

10.000 

10.0 

.042 

.6263 

Convergent Taper 

MPL 

.420 

3,880 

4.80 

,027 

.2557 


Table A. 4 Rotordynamic coefficients and leakage for the high-pressure 
turbine seal and the turbine interstage seal, 



































ORIGINAL PAGE [8 
OF POOR QUALITY 



K 

k 

C 

c 

lb /in 

lb /in 

lb sec/ln 

lb sec/ln 

Boost-Impeller 
Inlet (labyrinth) 

1206. 

740, 

.712 

4.43 X 10" 3 

Boost-Impeller 

Discharge 

(labyrinth) 

10,000. 

5,470. 

5.26 

63.3 X 10" 3 

Second-stage 
Turbine Tip Seal 

6,207. 

777 . 

.738 

22.3 X 10" 3 

First-stage 
Turbine Tip Seal 

5,739. 

568. 

.547 

20.3 X 10" 3 

Shrouded Inducer 
(one seal) 

62,000 

65,000 

98. 

0,0 


A. 5 Rocketdyne-calculated rotordynamic coefficients for labyrinth seals 
and shrouded-inducer selas at MPL; u = 20,218 cpra 



K 

•> 

C 

c 

lb /in 

lb /in 

lb sec/ln 

lb sec/in 

Boost-Impeller 
Inlet (labyrinth) 

2610. 

1760. 

1.08 

3 

5.79 X 10" 

Boost-Impeller 

Discharge 

(labyrinth) 

27,690. 

15,133. 

9.33 

131. X 10" 3 

Second-stage 
Turbine Tip Seal 

6,207. 

4 , 560 . 

2.93 

89.7 X 10" 3 

First-stage 
Turbine Tip Seal 

5,739. 

4,640. 

32.86 

92.9 X 10" 3 

Shrouded Inducer 
(one seal) 

154,000 

165,000 

166.0 

0.0 


A. 6 Rocketdyne-calculated rotordynamic coefficients for labyrinth seals 
and shrouded-inducer seals at FPL; U) = 30,381 cpm. 
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